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Abstract 

We present PAREVAL package containing a repository of theoretical physical mod- 
els used for (re-)evaluation of the fundamental physical constants (FPC). It holds all 
necessary data for building 105 (so called) observational equations and can be used 
in high precision calculations. Among repository models there are expressions for 
energy levels of hydrogen and deuterium (with 16 types of contributions), electron 
and muon magnetic moment anomalies, muonium ground-state hyperfine splitting, 
Zeeman energy levels in muonium. Each model is represented as Mathematica mod- 
ule with XML meta-data keeping information about the model (including data on 
dependence from others models). There are also modules for working with the basic 
FPC. 
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(i.e. expressions arising from physical theories) which can be written in terms of the 
basic fundamental physical constants (FPC). 

PAREVAL package was designed with following ideas in mind: (1) to simplify the 
procedure of the addition of new physical models ('breadth evolution'); (2) to sim- 
plify the procedure of the attachment of new models built from more general prin- 
ciples ('depth evolution'); (3) give an environment for comparison of the models. 
At present the models are precisely taken from [1,2] and the correctness of their 
representation in Mathematica was checked during our re-evaluation of the basic 
FPC. 
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[1] P. J. Mohr and B. N. Taylor, The 2002 CODATA Recommended Values of the 
Fundamental Physical Constants, Reviews of Modern Physics 77 1 (2005). 
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LONG WRITE-UP 



1 Introduction 

PAREVAL package is a set of Mathematica modules and notebooks which can 
be arranged by following: 

• The repository of physical models (i.e. Mathematica modules with some 
meta-data in XML); 

• Modules for the FPC usage. Contain data and functions for the FPC usage 
in high-precision calculations; 

• Notebooks with examples. The most interesting example is probably our 
re-evaluation of the basic FPC-2002; 

• Module for the final presentation of correlated quantities. Contain functions 
for the final rounding of results (according with techniques described in [3]). 

Each repository module defines some Mathematica symbols (variables or func- 
tions). Value of the defined symbol usually corresponds to some physical ex- 
pression. So the content of the repository we can represent as the realization 
of some phenomenological physical knowledge. 

Now the models are precisely taken from [1] and the correctness of their rep- 
resentation was checked during our re-evaluation of the basic FPC. 

Scientists from "Fundamental Constants Data Center" at National Institute 
of Standards and Technology (NIST) carried out a great work on collecting 
data which can be used for precision determination of the values of the FPC. 
We think that the theoretical expressions, on which the evaluation is based 
on, are interesting by itself and it's the main motivation for designing the 
PAREVAL package. 

New theoretical expressions can be easily add to the repository by the user. For 
example, one could find a module (Ml), which should be improved, introduce 
a new Mathematica symbol (or function) into it, place calculations for the 
symbol into separate module file and at last include the module name in 
dependence list in XML meta-data file of the first module Ml (thus it will be 
loaded by the corresponding repository function before the Ml module). We 
represent this to ourselves as 'breadth evolution' of the repository. 

Since Mathematica allows in principle to carry out any calculations, it's pos- 
sible to replace a model by another module calculated from more general 
principles. We call this by 'depth evolution'. 
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As one can see repository models can be written in terms of the basic funda- 
mental physical constants. Let's notice that the basic FPC is a some subset 
of the FPC which have following property: each of the FPC can be precisely 
expressed as the function of the basic ones. So after evaluation the basic FPC 
all the FPC can be calculated. 

1.1 A simple technology for accounting computational uncertainties 

In our FPC-2002 evaluation we have been generally based on the reviews 
[1], [2]. But in the realization of the models we embedded some our ideas. 
We've designed a simple technology for accounting for values with calculation 
uncertainties. 

Some parts of theoretical expressions for physical quantities come from nu- 
merical calculations, fits or approximations and only numerical values are 
usually known for them. They usually have uncertainties which should be ac- 
counted for. We solve the problem as following. Instead of use a scalar value 
A (which have an uncertainty U) we insert an object 'ss[A,U, Label]' (Mathe- 
matica function) into corresponding Mathematica expression and have created 
several function 'svalue', 'suncer' which suppress the 'ss[...]' objects (from the 
expression) and use corresponding information from it: 

• ' svalue [expr]' - gets out the value of the expression 'expr' in which 'A' value 
is used for every 'ss[A,U,Label]'; 

• 'suncer [expr]' - gets out the uncertainty of the expression 'expr' where ob- 
jects {'ss[A,U, Label]'} are considered as independent stochastic quantities 
(with corresponding normal distributions) with uncertainties {'U'}. 

'Label' data can be used to identify the object in the expression. 



2 The repository of physical models 

Each repository model is represented as Mathematica module (i.e. 'NAME.m' 
file, where 'NAME' is the name of the model) and 'NAME.xml' file with the 
following structure of XML tags: 

<m> <! — the model opening tag — > 

<dep> 'sequence of the names of models which should be loaded 

before this model ' </dep> 
<p> <! — parameters opening tag — > 

<inp> 'sequence of the input parameters ' </inp> 
<out> 'sequence of the output parameters' </out> 
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<sha> 'sequence of the shared parameters' </sha> 
</p> <! — parameters closing tag — > 
<title> 'the title of the model' </title> 
</m> <! — the model closing tag — > 

The repository models are listed in the Appendix A. 

These models were enough for building 105 (so called) observational equations 
(i.e. expressions which : (1) were measured experimentally (or in some special 
cases calculated or estimated theoretically); (2) can be expressed in terms of 
the FPC) selected by the CODATA-2002. 

These observational equations can be separated into two groups "A" (princi- 
ple observational equations for determination of the Rydberg constant) and 
"B" (principle observational equations for determination of some others basic 
FPC). 

In order to show the correctness of the repository models we prepared several 
web-pages with the data on the observational equations 1 . 

Firstly, we would like to present a table where experimental data are compared 
with theoretical ones: 

• Data on the observational equations. Group "A" http://hermes.ihep.su/fpc/table-oes-a.xhtml; 

• Data on the observational equations. Group "B" http:/ /hermes. ihep.su/fpc/table-oes-b.xhtml; 

Let's make some comment: 

• Correlation coefficients of the input data for FPC-2002 evaluation are pre- 
sented on the page: http:/ /hermes. ihep.su/fpc/table-oes-corr.xhtml 

• "ID" column: The numbers after A'-, B'- indicate on the position from the 
beginning of the group (For "A" group it corresponds to CODATA-2002 
"LSA index" (see, for example, a table from NIST web site 
http://physics.nist.gov/cuu/LSAData/qdata.pdf), "B'l" corresponds to 51 
"LSA index", "B'2" to 52 and so on). So ordering of the observational 
equations is precisely the same as used in the CODATA-2002 evaluation); 

• "Exp. uncer.", "Rel. exp. una": experimental and relative experimental 
uncertainties correspondingly; 

• "Theor. uncer.", "Rel. theor. una": theoretical and relative theoretical un- 
certainties correspondingly. Theoretical uncertainty arises from the uncer- 
tainties of the FPC (including (^-corrections due to uncalculated terms) ac- 
cording to the law of uncertainty propagation; 



Web-browser should support MathML format to display the 
pages. For Internet Explorer you can download MathPlayer plug-in: 
http : / / www . dessci .com/en/ products / mat hplayer / . 
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• "Calc. uncer.": Calculational uncertainty was shortly discussed in the sec- 
tion 1.1. It's due to the uncertainties that some of auxiliary parameters in 
theoretical expressions have. Large calculation uncertainty for A'l observa- 
tional equation is mainly due to uncertainty of C50 coefficient from expres- 
sion for three photon correction ([1], p. 82). We believe that it was reduced 
by introduction the <5n(h, L,j) 2 with corresponding auxilary observational 
equations (and corresponding correlation coefficients). From other hand it 
shows that computational uncertainty is rough estimation of the uncertainty 
(but possibly more reliable); 

• "Cal.uncer/Th. uncer" : Ratio of the computational uncertainty to the the- 
oretical one. It's is an indicator: usually it should be less then 0.1; 

• " (Exp.-Theor.) /Una" is an indicator of agreement of theoretical and experi- 



mental values. "Una" is the total uncertainty: Una = Ju^ xperiment + Uf heory 

(where Ux is corresponding uncertainty); 

• " Original expression" is original expression for observational equation from 
the repository model. 

Secondly, we show tables of data calculated by NIST (after FPC-2002 evalu- 
ation) and with the PAREVAL: 

• Comparison values of the observational equations of the group "A" : CODATA- 
2002/NIST vs. PAREVAL 

http: / /hermes. ihep.su/fpc/table-oes-a_NIST_PARE.xhtml; 

• Comparison values of the observational equations of the group " B" : CODATA- 
2002/NIST vs. PAREVAL 

http://hermes.ihep.su/fpc/table-oes-b_NIST_PARE.xhtml. 
Comments are followings: 

• Rows are sorted by the values of fifth column - indicator of computational 
agreement; 

• "Exp. Uncer." is experimental uncertainty. Usually it defines the level of 
required precision so we compare the values via it: it should be less then 10% 
normally (in our case the maximum value is less then 7% with maximum 
for AT4 item). 

So as one can see we have good agreement with NIST calculations of the 
observational equations. 



2 corrections to the energy levels of the hydrogen and deuterium due to uncalculated 
terms 
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3 The results of our FPC-2002 re-evaluation 



The most interesting example of the repository usage is probably our re- 
evaluation of the basic FPC-2002. 

We obtained good agreement with CODATA-2002 values of the basic FPC: 
the maximum value of the 'Pull' = \{z PA reval - zcodata)/U PA reval\ is less 
then 4% (where by zx we denote corresponding value of the basic FPC). 

All the constants uncertainties obtained by us are very closely to CODATA- 
2002 ones in the meaning | U pareval / U cod at a — 1| < 1% , where by Ux we 
denoted uncertainty of the constant (for CODATA-2002 and our results). 

One can see results of our evaluation from the followings web-pages 3 : 

• The table with the basic FPC http://hermes.ihep.su/fpc/table-fpc.xhtml; 

• The correlation matrix of the basic FPC http: / /hermes. ihep.su/fpc/table-fpc-corr .xhtml; 

To determine the number of significant figures for final rounding procedure 
for the uncertainties of the basic FPC we use an expression [3]: 



, where A^ in is the minimal eigenvalue of the correlation matrix of the basic 
FPC. So we got P^ 1 = 4 (as one can see from the table of the basic FPC from 
the web-page). 

In order to save correlation matrix to be positive defined after final rounding 
procedure we use the following expression for number of significant figures: 



We got = 7. 



4 Conclusion 

By means of PAREVAL package the modern theoretical expressions which 
were collected by 'Fundamental Constants Data Center" at NIST become 

3 Web-browser should support MathML format to display the 
pages. For Internet Explorer you can download MathPlayer plug-in: 
http : / / www . dessci .com/en/ products / mat hplayer / . 




[3]: 
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accessible for scientific community as repository of physical models in Mathe- 
matica computer algebra system. So the repository represents a reference-book 
of standard physical expressions conjugated with calculation media. This al- 
lows to see, calculate and compare theoretical expressions with simplicity. 
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6 Appendix A 



The repository modules are listed in following two tables 4 : 



4 All indicated pages in the tables refer usually to [1] if it's not written anything 
else. 



Title 


Output Symbol (s) 


Module file 


Electron magnetic moment anomaly 
(p.474-476) 


a» e » [a, dae] 


a-e.m 


Muon magnetic moment anomaly 
(p. 86-89) 


a>y [a, damn] 


a-mu.m 


Models for energy levels of hydrogen and deuterium atoms: 


- General contributions (p. 77-84) 


ElbyParts Eltot 


e-l-2.m 


- Two-photon corrections (p. 80) 


EU 


e-L4.m 


- Three-photon corrections (p. 82) 


EU 


e-L6.m 


- Finite nuclear size contribution 
(p.82) 


EIns 


e-Lfns.m 


- Relativistic recoil non-leading 
terms contribution (p. 78) 


tL 1 F?ni 


e-l_rr2.m 


- Radiative-recoil correction 
(p.83) 


EIrr 


e-l_rr.m 


Ratios of bound-particle to 
free-particle g-f actors (p. 93) 


rgeHge rgpHgp 
geDge gmupMugmup 
geMuge gdDgd 


g-fact-rat.m 


Bound-state g-factors (for bound 
electron in 12 C 5+ and 16 7+ ) 
(p.89) 


gg"C" [5] 
gg-O" [$] 
9»e» 


g-fact.m 
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Title 


Output Symbols 


Module file 


Muoniom ground-state hyperfine 
splitting (p.94-96) 


dnMUth 


mu-hs.m 


Zeeman energy levels in muonium 
([2],p.386-387; [l],p.94-96) 


v[fp\ 


mu-zee.m 


Ground-state ionization energies 
(p. 11, table IV) 


matrix4 


table-e_b.m 


Set of the observational equations 
for data "A" for the FPC-2002 
evaluation (p. 57) 


OEListDefAl 
OEListDefDelta 


tdata-a.m 


Set of the observational equations 
for data "B" for the FPC-2002 
evaluation (p. 59-60) 


OElistDefB OElist 


tdata-b.m 


Mathematica variables for 
calculations with the FPC-2002 
(the results of our re-evaluation) 


varsFPC subsFPC 
uncerFPC unitsFPC 
namesFPC 
corrFPC covFPC 
texsymbFPC 


reFPC-2002a.m 


Functions for accounting for 
computational uncertainties 


svalue[expr] 
suncer[expr] 


nnum.m 
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